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Abstract. We continue the study of the free energy of quantum lattice spin 
systems where to the local Hamiltonian H an arbitrary mean field term is 
added, a polynomial function of the arithmetic mean of some local observables 
X and Y that do not necessarily commute. By slightly extending a recent 
paper by Hiai, Mosonyi, Ohno and Petz [3], we prove in general that the free 
energy is given by a variational principle over the range of the operators X and 
Y. As in [9], the result is a noncommutative extension of the Laplace- Varadhan 
asymptotic formula. 



1. Introduction 

1.1. Large deviations. One of the highlights in the combination of analysis and 
probability theory is the asymptotic evaluation of certain integrals. We have here 
in mind integrals of the form, for some real- valued function G, 

dfj, n (x) exp{v n G(x)}, u„ / +00 as n y +00 (1-1) 

for which the measures fi n satisfy a law of large numbers. Such integrals can 
be evaluated depending on the asymptotics of the The latter is the subject 
of the theory of large deviations, characterizing the rate of convergence in the 
law of large numbers. In a typical scenario, the /i n are the probabilities of some 
macroscopic variable, such as the average magnetization or the particle density in 
ever growing volumes v n and as distributed in a given equilibrium Gibbs ensemble. 
Then, depending on the case, thermodynamic potentials J make the rate function 
d/j, n (x) ~ da;exp{— v n J{x)} in the sense of large deviations for Gibbs measures, 
see [T^J [H [HJ [52]. That theory of large deviations is however broader than 
the applications in equilibrium statistical mechanics. Essentially, when the rate 
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function for fj, n is given by 3 ', then the integral (jl.ip is computed as 

— log [ dn n (x)exp{v n G(x)} — > sup{G(x) - J(x)} (1.2) 

V n J nS+oo x 

This is a typical application of Laplace's asymptotic formula for the evaluation of 
real-valued integrals. The systematic combination with the theory of large devia- 
tions gives the so called Laplace- Varadhan integral lemma. 

We first recall the large deviation principle (LDP). Let (M, d) be some complete 
separable metric space. 

Definition 1.1. The sequence of measures [i n on M satisfies a LDP with rate 
function J : M — ► R+ U {+00} and speed v n E K+ if 

1) J is convex and has closed level sets, i.e., 

{J~\x), x<c} (1.3) 

is closed in (M, d) for all c <E R+; 

2) for all Borel sets U C M with interior int U and closure cl U, one has 

liminf — log// n (t/) > — inf J{u) 

n/"'+oc V n u£int U 

limsup — \og[i n (U) < — inf J{u) 

We say that the rate function J is good whenever the level sets (jl.3p are compact. 

For the transfer of LDP, one considers a pair (ji n , v n ), n /* 00 of sequences of 
absolutely continuous measures on (Af, d) such that 

— —(x) = exp{w„ G(x)}, fi n — almost everywhere 

for some measurable mapping G : M — > M. We now state an instance of the 
Laplace- Varadhan lemma. 

Lemma 1.1 (Laplace- Varadhan integral lemma). Assume that G is bounded and 
continuous and that the sequence (fJ> n ) satisfies a large deviation principle with good 
rate function J and speed v n . Then (u n ) satisfies a large deviation principle with 
good rate function G — J and speed v n ■ 

For more general versions and proofs we refer to the literature, see e.g. [2TJ [22l 
H3E1I2]; it remains an important subject of analytic probability theory to extend 
the validity of the variational formulation (11.2[) and to deal with its applications. 

1.2. Mean-field interactions. From the point of view of equilibrium statistical 
mechanics, one can also think of the formula (|l.ip as giving (the exponential of) 
the pressure or free energy when adding a mean field type term to a Hamiltonian 
which is a sum of local interactions. 

The choice of the function G is then typically monomial with a power decided 
by the number of particles or spins that are in direct interaction. For example, the 
free energy of an Ising-like model with such an extra mean field interaction would 
be given by the limit 

}%m l °z E -p(-/^aW + a p |a|(^$>) p ) (1.4) 

1 1 ve{+, -} A 1 1 ieA 

for p — 1,2,..., where H^(rf) is the (local) energy of the spin configuration 77 
and the limit takes a sequence of regularly expanding volumes A to cover some 
given lattice. The case p = 1 corresponds to the addition of a magnetic field 



3 



Ai; p = 2 is most standard and adds effectively a very small but long range two- 
spin interaction. Higher p— values are also not uncommon in the study of Ising 
interactions on hypergraphs, and even very large p has been found relevant e.g. in 
models of spin glasses and in information theory [2] . 
The form (|l.ip is easily recognized in (|1.4|) . with 

Hn(x) ~ ^2 exp{-f3H A (r])}, v n = |A| 

*?e{+ ,-} A ,E ieA >7«=s|A| 

and the function G(x) = X p x p . The Laplace- Varadhan lemma applies to (|1.4j) since 
we know that the sequence of Gibbs states with density ~ exp{— (3Ha(-)} satisfies 
a LDP with a good rate function J7" c i and speed |A|. The result reads that (|1.4|) is 
given by the variational formula 

sup {\ p uP ~ J cl (u)} (1.5) 
ue[-i,i] 

In noncommutative versions the local Hamiltonian H and the additional mean 
field term are allowed not to commute with each other. That is natural within the 
statistical mechanics of quantum spin systems and this is also the context of the 
present paper. 

1.3. Noncommutative extensions. Although it has proven very useful to think 
of integrals (jl.ip within the framework of probability and large deviation theory, it is 
fundamentally a problem of analysis. However, without such a probabilistic context, 
the question of a noncommutative extension of the Laplace- Varadhan Lemma 11.11 
becomes ambiguous and it in fact allows for different formulations, each possibly 
having a physical interpretation on its own. 

One approach is to ask for the asymptotic evaluation of the expectations 

lim ^-logw A (e |A|G( ^ A) ) (1.6) 
Ayz d |A| v 7 

under a family of quantum states uja where X\ would now be the arithmetic 
mean of some quantum observable in volume A. To be specific, one can take 
wa ~ exp{— f3H\} a quantum Gibbs state for a quantum Hamiltonian H\ and 
X\ = {^2 ie \Xi)/\A\ the mean magnetization in some fixed direction. Arguably, 
this formulation is closely related to the asymptotic statistics of outcomes in von 
Neumann measurements of X\, [13] . \12 \ fl4| . IPO] . 

A more general class of possible extensions is obtained by considering the limits 

of 

-^logTr^Je^^A))*, A/Z rf (1.7) 

for different K > 0, where <ta is the density matrix of a quantum state in volume A. 
For the canonical form a\ = exp(—f3H\)/Z^ with local Hamiltonian Ha at inverse 
temperature 0, (| 1 . 7|) becomes 

-L log * Tr A (e-*** , A / Z" (1.8) 

There is no a priori reason to exclude any particular value of K from consideration. 
Two standard options are: K = 1, which corresponds to the expression (|1.6[) above, 
and K /* +oo, which, by the Trotter product formula, boils down to 

i- log ^ Tr A (e-^+l A l G <^>) , A / Z d (1.9) 

I I Z A 

which is the free energy of a corresponding quantum spin model, cf. (|1.4[) . In the 
present paper, we study the case K /* +oo (without touching the question of 
interchangeability of both limits). 
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One of our results, Theorem 13.11 with Y = Y A = 0, is of the form 

lim J-logTrAfe-^+l^ ^) = sup {G(u) - J(u)} (1.10) 

A/Z d |A| -||X||<«<||X|| 

Note that we omitted the normalization factor 1 /Z^ since it merely adds a constant 
(independent of G) to (|1.9p . In the usual context of the theory of large deviations 
formula (|1.10[) arises as a change of rate function. However, while our result (TlXO)) 
very much looks like Varadhan's formula in Lemma ll.li there is a big difference in 
interpretation: The function J is not as such the rate function of large deviations 
for X A . Instead, it is given as the Legendre transform 

J(u) = sup {tu - q(t)}, ueR (1.11) 

of a function q(-) which is the pressure corresponding to a linearized interaction, 
i.e. 

q(t) = lim -i-logTr A (e-^+*l A l^)) (1.12) 

A/TL d |A| v ' 

1.4. Several non-commuting observables: Towards joint large deviations? 

In the previous Section fl.3[ we made the tacit assumption that there is a single ob- 
servable X A corresponding to some operator on Hilbert space. However, in formula 
11.41 the observable m J2n£ A Vi could equally well represent a vector- valued mag- 
netization which, upon quantization, would correspond to several non-commuting 
observables X A ,Y A , say, the magnetization along the x and y-axis, respectively In 
the commutative theory, this case does not require special attention; the framework 
of large deviations applies equally regardless of whether the observable takes values 
in R or R 2 . Obviously, this is not true in the noncommutative setting and in fact, 
we do not even know a natural analogue of the generating function (|1.6p . since 
we do not dispose of a simultaneous Von Neumann measurement of X A and Ya. 
One can take the point of view that this is inevitable in quantum mechanics, and 
insisting is pointless. Yet, as A /* Z d , the commutator 

[X A ,Y A } = 0(-^-) (1.13) 

vanishes and hence the joint measurability of X A , Ya is restored on the macroscopic 
scale. We refer the reader to [Tg] where this issue is discussed and studied in more 
depth. 

The advantage of the approach via the Laplace- Varadhan Lemma is that one can 
set aside these conceptual questions and study 'joint large deviations' of X A and 
Ya by choosing G to be a joint function of X A and Ya, for example a symmetrized 
monomial 

G(X A ,Y A ) = (X A ) k (Y A ) l + (X A ) l (Y A ) k , for some k,leN (1.14) 

and check whether the formula ()1.10p remains valid with some obvious adjustments. 
This turns out to be the case and it is our main result: Theorem 13. II 

1.5. Comparison with previous results. The asymptotics of the expression ()1.9|) 
was first studied and the result (| 1 . 10[) was first obtained by Petz et al. [M], in the 
case where the Hamiltonian H A is made solely from a one-body interaction. The 
corresponding equilibrium state is then a product state. In [5], Hiai et al. gener- 
alized this result to the case of locally interacting spins but the lattice dimension 
was restricted to d — 1. However, the authors of [9] argue that the restriction to 
d = 1 can be lifted in the high-temperature regime. The main reason is that their 
work relies heavily on an asymptotic decoupling condition which is proven in that 
regime, 1 . One should observe here that this asymptotic decoupling condition in 



fact implies a large deviation principle for X\, as follows from the work of Pfistcr 
[T7] , Hence, in the language of Section fT~3l 9 evaluates (|1.9|) (the case K = oo) in 
those regimes where the (|1.6[) (the case K = 1) is analytic. 

The present paper elaborates on the result of [9] in two ways. First, we remark 
that, in our setup, the decoupling condition is actually not necessary for (|1.10j) to 
hold, and therefore one can do away with the restriction to d = 1 or high tempera- 
ture. Hence, again referring to Section fl. 31 the case K = oo can be controlled even 
when we know little about the case K = 1. To drop the decoupling condition, it 
is absolutely essential that we start from finite-volume Gibbs states, and not from 
finite- volume restrictions of infinite- volume Gibbs states, as is done in [9]. 

Second, we show that by the same formalism, one can treat the case of several 
noncommuting observables, as explained in Section 11.41 The most serious step in 
this generalization is actually an extension of the result of |16j to noncommuting 
observables. This extension is stated in Lemma 16. II and proven in Section [3 

Note While we were finishing this paper, we learnt of a similar project by W. de 
Siqueira Pedra and J-B. Bru. Their result [3] is nothing less than a full-fledged 
theory of equilibrium states with mean-field terms in the Hamiltonian, describing 
not only the mean-field free energy (as we do here), but also the states themselves. 
Also, their results hold for fermions, while ours are restricted to spin systems, and 
they provide interesting examples. Yet, the focus of our paper differs from theirs 
and our main result is not contained in their paper. 

1.6. Outline. In Section[2j we sketch the setup. We introduce spin systems on the 
lattice, noncommutative polynomials and ergodic states. Section [3] describes the 
result of the paper. The remaining Sections @1 [SJ [B] and [7] contain the proofs. 



2.1. Hamiltonian and observables. We consider a quantum spin system on the 

regular lattice Z d , d = 1, 2, We briefly introduce the essential setup below, and 

we refer to |lll 119] for more expanded, standard introductions. 

The single site Hilbert space TL is finite-dimensional (isomorphic to C n ) and for 
any finite volume A C Z d , we set Ha = ®aH. The C*-algebra of bounded operators 
on Ha is denoted by Ba = B(Ha)- The standard embedding Ba C Ba> for Ac A' 
is assumed throughout. The quasi- local algebra il is defined as the norm closure of 
the finite-volume algebras 



Denote by t^, i G Z d , the translation which shifts all observables over a lattice 
vector i, i.e., 73 is a homomorphism from Ba onto Bi+A- 

We introduce an interaction potential $, that is a collection ($a) of Hermitian 
elements of Ba, labeled by finite subsets icZf We assume translation invariance 
(i) and a finite range (ii): 

i) r,($ A ) = $ l+j4 for all finite A C 7L d \ 

ii) there is a d max < oo such that, if diam(A) > d max , then <&a = 0. 
In estimates, we will frequently use the number 



2. Setup 




(2-1) 



A finite 




(2.2) 



A30 



The local Hamiltonian in a finite volume A is 




(2.3) 



A<ZA 
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which corresponds to free or open boundary conditions. Boundary conditions will 
however turn out to be irrelevant for our results. We will drop the superscript $ 
since we will keep the interaction potential fixed. 

Let X,Y, . . . denote local observables on the lattice, located at the origin, i.e., 
SuppX (which is defined as the smallest set A such that X g Ba) is a finite set 
which includes € Z d . 

We write 

X A := T > X ( 2 - 4 ) 

jez d ,Supp-rjXcA 

and 

x a •= jx\ Xa ( 2 - 5 ) 

for the corresponding intensive observable (the 'empirical average' of X). 

All of these operators are naturally embedded into the quasi-local algebra il. At 
some point, we will also require the infinite volume intensive observable 

X ~ Xa/oo 

Some care is required in dealing with X since it does not belong to the quasi-local 
algebra 11. We will further comment on this in Section 12.31 

2.2. Noncommutative polynomials. We will perturb the Hamiltonian ff * by a 
mean field term of the form |A|G(Xa, Y a ) where G is a 'noncommutative polyno- 
mial' of the operators Xa, Ya, e.g. as in (|1.14[l . 

In this section, we introduce these noncommutative polynomials G as quantiza- 
tions of polynomial functions g. First, we define 

Ran(X,Y):=[-||X||,||X||]x[-||y||,||y||] (2.6) 

This definition is motivated by the fact that ('sp' stands for spectrum) 

spX A x spY A C Ran(X, Y), for all A (2.7) 

Let g be a real polynomial function on the rectangular set Ran(X, Y). Using the 
symbol I for the collection of all finite sequences from the binary set {1,2}, any 
map G : X — > C is called a quantization of g whenever 

N 

^2 G{a)x a(1) ...x a{n) =g(x 1 ,x 2 ) (2.8) 

n=0 a=(a(l),...,Q(n))£l 

for all (xi,X2) G Kan(X, Y) and for some N G N. A quantization G is called 
symmetric whenever 



G(a(l), . . . , a(n)) - G(a(n), . . . , a(l)). (2.9) 
Any such symmetric quantization G defines a self-adjoint operator 

N 

G(X,Y) = J2 G{o)X a{1) ...X a(n) (2.10) 

n=0 a=(a(l),...,a(n))el 

taking X 1 = X and X 2 = Y. 

In the thermodynamic limit, one expects different quantizations of g to be equiv- 
alent: 

Lemma 2.1. Let G andG' be any two quantizations of g : Ran(A, Y) — > K. Then 
\\G(X A ,Y A )-G'(X A ,Y A )\\ < C( *' Y) (2.11) 
for some C(X, Y) < oo, and for all finite volumes A. 
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Proof. This is a simple consequence of the fact that the commutator of macroscopic 
observables vanishes in the thermodynamic limit, more precisely, 

IPa,?a]|| < Su PP X| x ||Y|||Sup P y|. (2.12) 

□ 

Indeed, our results, Theorems 13.11 and 13.21 do not depend on the choice of quan- 
tization. This can also be checked a priori using the above lemma and the log-trace 
inequality in (|3.10p . 

2.3. Infinite-volume states. A state wa is a positive linear functional on B\, 
normalized by \\uj\\\ = = 1. An example is the tracial state, wa(-) ~ Tr A (-)- 

In general we consider states lo\ as characterized by their density matrix a\, wa(-) = 
Tr A (<7 A -)- 

An infinite volume state wisaa positive normalized function on the C*-algebra 
il (the quasi-local algebra). It is translation invariant when oj(A) = uifaA) for all 
j G Z d and A € il. A translation-invariant state ui is ergodic whenever it is an 
extremal point in the convex set of translation invariant states. A state is called 
symmetric whenever it is invariant under a permutation of the lattice sites, that is, 
for any sequence of one-site observable A\, . . . ,A n G £>{ } C il and i\, . . . , i n G Z d 

Lu{ Til {Ax)T i2 (A 2 ) ...r in {A n )) = u;(r iir(1) (A^n^ (A 2 ) . . . r iw(n) (A n )) (2.13) 

for any permutation 7r of the set {1, . . . ,rt}. The set of ergodic, resp. symmetric 
states on il is denoted by «Serg,S sym , respectively. 

At some point we will need the theorem by St0rmer [20] that states that any 
lo G <S sym can be decomposed as 

&„{<t>)<t> (2.14) 

prod. 

for some probability measure u u on the set of product states. Of course, the set of 
product states can be identified with the (finite-dimensional) set of states on the 
one-site algebra B^ j = B(TL). 

For a finite- volume state o->a on Ba, we consider the entropy functional 

S'(wa) = Sa(wa) = -TroAlogo-A (2-15) 
The mean entropy of translation-invariant infinite-volume state u> is denoted as 

s(uj) :— lim -r-—rS(ui\), with c<ja ^1 K (restriction to A) (2.16) 
Ayz d |A| IBa 

In this formula and in the rest of the paper, the limit lim A ^ Z d is meant in the sense 
of Van Hove, see e.g. [TT] |TS] . Standard properties of the functional s are its affinity 
and upper semicontinuity (with respect to the weak topology on states). 

In Section |2~T| we mentioned the 'observables at infinity' X and Y, postponing 
their definition to the present section. Expressions like u(X l Y k ) (for some positive 
numbers I, k) can be defined as 

co(X l Y k ) := lim (2.17) 

A,A'/"Z d 

provided that the limit exists. We use the following standard result that can be 
viewed as a noncommutative law of large numbers 

Lemma 2.2. For lo G <S crg , the limit (|2.17p exists and 

iu(X l Y k ) = [uj(X)] l [uj(Y)] k (2.18) 
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Note that u>(X) = u>(X) and u)(Y) — lu{Y) by translation invariance. An imme- 
diate corollary is that for a noncommutative polynomial G which is a quantization 
of g (see Section |2~2|) . and ui e <S crg : 

u{G{X,Y))=g(w(X),w{Y)) (2.19) 

For the convenience of the reader, we sketch the proof of Lemma l2.2l in Appendix 

El 

Finally, we note that Lemma 12.21 does not require the state to to be trivial 
at infinity. Triviality at infinity is a stronger notion which is not used in the 
present paper. In particular, the state p, constructed in Section U] is ergodic, but 
not trivial at infinity, since it fails to be ergodic with respect to a subgroup of lattice 
translations. 

3. Result 

Choose X, Y to be local operators and let i7* be the Hamiltonian corresponding 
to a finite-range, translation invariant interaction as in Section ^. II Let G(-, •) be 
a symmetric quantization of a polynomial g on the rectangle Ran(X, Y), as defined 
in Section [2~2l and define the "G-mean field partition function" 

Z«($) := TT A (e- H ^ G ^^) (3.1) 

with X\, Y\ empirical averages of X, Y . The following theorem is our main result: 

Theorem 3.1. Define the pressure 

p(u,v)= lim -5-logTr A e- H A+^A+^A ( 32) 

A/"Z d |A| 

and its Legendre transform 

I(x,y) — sup (ux + vy — p(u, v)) (3-3) 

0,u)eR 2 

Then 

lim m logZf($)= sup (g(x,y) - I(x,y)) (3.4) 

Ayz d |A| (x,y)m 2 



where the limit A / Z is in the sense of Van Hove, as in (13. 2|) . In particular, the 
LHS of (|3.4p does not depend on the particular form of quantization taken. 



As discussed in Section [IJ our result expresses the pressure of the mean field 
Hamiltonian through a variational principle. To derive this result, it is helpful to 
represent this pressure first as a variational problem on a larger space, namely that 
of ergodic states, as in Theorem 13.21 Theorem 13.11 follows then by parametrizing 
these states by their values on X and Y. 

We also need the 'local energy operator' associated to the interaction $ as 

^ : =Enj^- ( 3 - 5 ) 

Theorem 3.2 (Mean- field variational principle). Let s(-) be the mean entropy func- 
tional, as in Section \2.!A Then 

lim JL bg ($) - sup (g(u{X), u(Y)) + a(w) - u{E 9 )) (3.6) 
A/^ |A| w6Serg 



To understand how the first term on the RHS of (|3.6| originates from (|3.1|) . we 
recall the equality (|2.19p for ergodic states lo. 

The proof of Theorem l3.2l is postponed to Sections [5] and [S] Here we prove that 
Theorem 13.11 is a rather immediate consequence of Theorem [ 
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Proof of Theorem \3.1\ Writing the right-hand side of (|3.6[) in the form 



sup (g(x,y) - I(x,y)) (3.7) 

0,y)GR 2 



where 



I(x,y)= inf (-s(u)+w(£$)) (3.8) 

w(X)=<o, w(V)=y 

is a convex function on R 2 , infinite on the complement of Ran(A, Y). It is lower 
semi-continuous by the lower semi-continuity of — s. By using the infinite- volume 
Gibbs variational principle [HI [11] , its Legendre-Fenchel transform reads 

sup (ux + vy — I(x,y)) = sup (s(u>) — cj(.E$) + uu)(X) + vu)(Y)) 
(x,i/)eR 2 wes mg 

= p(u,v) (3.9) 

The equality 1 = 1 then follows by the involution property of the Legendre-Fenchel 
transform on the set of convex lower-semicontinuous functions, see e.g. [19] . □ 



Independence of boundary conditions. Observe that both Theorem 13.11 
and Theorem [321 have been formulated for the finite volume Gibbs states with open 
boundary conditions. It is however easy to check that this choice is not essential and 
other equivalent formulations can be obtained. Indeed, by the standard log-trace 
inequality, 

lo g Tr A (e-^ A + WA +l A l G ^ A ^ A )) -logTr A (e-^ A+ l A l G ^ A ' ?A >) < ||W A || (3.10) 

and hence if one chooses W\ such that lim A ^ Z d ||Wa||/|A| = 0, then we can replace 
-0H\ by -f3H A + W\ in Theorem O and Theorem EU 

One also sees that it suffices to prove Theorem 13.11 and Theorem 13. 21 for the case 
that A = (a[— L,L]) d and a /* oo. Convergence for A /* 1 d in the sense of Van 
Hove then follows by the above log-trace inequality. 

Finite-range restrictions. It is obvious that our paper contains some restric- 
tions that are not essential. In particular, by standard estimates (in particular, 
those used to prove the existence of the pressure, see e.g. [T!5]) one can relax the 
finite-range conditions on the interaction <E> to the condition that 

eW<-. (3-11) 

.430 1 1 

and similarly for the local observables X, Y. Moreover, it is not necessary that G 
is a noncommutative polynomial. Starting from (|3.10|) . one checks that it suffices 
that G can be approximated in operator norm by noncommutative polynomials. 



4. Approximation by ergodic states 

In this section, we describe a construction that is the main ingredient of our 
proofs, as well as of those in [3] and [TB]. This construction will be used in Sections 
[SandE 

Let V be a hypercube centered at the origin, i.e., V = [-L, L] d for some L > 1 
and let 

dV := {i G V | 3i' 6 Z d \ V such that i, i' are nearest neighbours} (4-1) 
Consider a state \iy on By- 

We aim to build an infinite- volume ergodic state out of \iy . First, we cover the 
infinite lattice Z d with translates of the hypercube V such that II 1 = U ie ( iZ )dV + i. 
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Then, we define the block product state 

p := ® [iv ° n (4.2) 

i<E(ZL) d 

where \iy o r, is a state (translate of /iy) on Si+y. We define also the 'translation- 
average' of p, 

1 1 jev 

We can now check the following properties: 

We have the exact equality of entropies 

s{p) = m = TtfS(»v) ( 4 - 4 ) 

This follows from the affinity of the entropy in infinite-volume. A remark 
is in order: A priori, the infinite-volume entropy is defined for translation- 
invariant states, whereas p is only periodic. However, one easily sees that 
the entropy can still be defined, e.g. be viewing p as translation-invariant 
with a bigger one-site space, and correcting the definition by the 'volume' 
of this one-site space. 

The state p is ergodic. This follows, for example, from Proposition 1.7.9 in 
[T9"] . as is also shown in [5] via an explicit calculation. Note however that 
the p is in general not ergodic with respect to the translations over the 
sublatticc (LZ) d . This phenomenon (though in a slightly different setting) 
is commented upon in [19] (the end of Section III. 5). 

The state p is a good approximation of fj,y for observables which are em- 
pirical averages, provided V is large. Consider the local observable X 
as in Section 12.11 A translate TjX can lie inside a translate of V, i.e. 
SuppTjX C V + i for some i £ X d /V = (LZ) d , or it can lie on the bound- 
ary between two translates of V. The difference between p(X) = p(X) and 
fiv(Xv) clearly stems from those translates where X lies on a boundary, 
and the fraction of such translates is bounded by 

|SuppX|x^ (4.5) 

Hence 

- , \dV\ 
\p(X) - vv(X v )\ < ||X|[|SuppX| x 1— 1 (4.6) 



\V\ 



and also \p(X) - p{X A )\ < \\X\\\ SuppX| x J^l. 



5. The lower bound 
In this section, we prove the following lower bound. 
Lemma 5.1. Recall Z^(^) as defined in (13. ip . Then 

liminf-|-log^($) > sup ((g(w(X),u;(Y))+8(u)-u(E*)) (5.1) 

where all symbols have the same meaning as in Section^ 
Proof. Consider a state ui € S CYg . We show that 

liminf ^- log Z%(<1>) > g(Lu(X),Lu(Y)) + s(lu) - lu(Ez) (5.2) 
Ayz d |A| 
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Consider, for each volume A, the restriction uj a := u>\ B . By the finite- volume 
variational principle 

±\ogZ°(<P) > u> A (G(X A ,Y A )) + ynS(cj A ) - ^a^a) (5.3) 

The following convergence properties apply (A f Z d in the sense of Van Hove): 
1) 

w a (G(X a ,Y a ))=lo(G(X a ,Y a ))^ 9 (lj(X),lj(Y)), A / Z d (5.4) 

2) 

j^jS(wA) - s(u>), k/7L d (5.5) 

3) 

^ u>{H A ) -» u(E 9 ), A/Z d (5.6) 

The relation (|5.6|) is obvious from the summability of <£>, see Section 12.11 The 
convergence ()5.5[) is the definition of the mean entropy s. Finally, (|5.4j) follows 
from the ergodicity of u) as explained in Section 12.31 

The relation (|5.2p now follows immediately, since one can repeat the above con- 
struction for any ergodic state uj. □ 

6. The upper bound 

6.1. Reduction to product states. In this section, we outline how to approxi- 
mate 

pqlogZf(S) (6.1) 

by a similar expression involving the partition function of a block-product state. 
Fix a hypercube V = [-L, L] d and cover the lattice with its translates, as explained 
in Section [3J From now on, A is chosen such that it is a 'multiple' of V, which is 
sufficient by the remark at the end of Section [3] Define the observables 

HX=H*> V , XX, YX (6.2) 

by cutting all terms that connect any two translates of V, i.e., 

3i 6 1 d /V : Supp TjAcV + i 

and analogously for HX and Y A - One can say that these observables with super- 
script V are 'one-block' observables with the blocks being translates of V. One 
easily derives that 

11*1 - X A \\ < \\X\\\ SuppX|H, \\HX -H A \\< r($)|A|ffl (6.4) 

with the number r($) as defined in Section [2.11 
Using the log-trace inequality, we bound 

^logTr A ( e -^A+|A|G(x A ,? A )) _ i_ logXrA (e-BX+\MG(xX,YZ^ (6 . 5) 
as follows 

m < ^\\h a -hX\\ + \\g(X Ai Y a )-g(XX,yX)\\ 
< (r($) + c ff (||x|||Su PP x| + ||y|||Su PP y|))H 
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where C g is constant depending on the function G. The second term of (|6.5p is 
clearly the pressure of a product state with mean field interaction. We will find an 
upper bound for this pressure by slightly extending the treatment of Petz et al. in 
[l6j . We prove an 'extended PRV -lemma, Lemma l6~Tl in the next section. 

6.2. The extended Petz-Raggio-Verbeure upper bound. In this section, we 
outline the bound from above on the quantity 

logTr A (e-Hl+WG(xX,Y^ (6 6) 

that appeared in (|6.5|) . 

To do this, let us make the setting slightly more abstract. Consider again the 
lattice % d with the one-site Hilbert space Q, which should be thought of as 

Q := ® V H (6.7) 

In words, the sites of the new lattice are actually blocks of the old lattice. In 
particular, the trace TrA has a slightly different meaning as before since the one- 
site Hilbert space has changed, i.e. the trace is on the Hilbert space ®aQ instead 
of ®hH- Let D,A,B be one-site observable on the new lattice, i.e. D,A,B are 
Hermitian operators on Q. The extended PRV (Petz-Raggio-Verbeure) states that 

Lemma 6.1 (Extended PRV). Let all symbols have the same meaning as in Sec- 
tions \KJS2^ME3i except that the one-site Hilbert space is changed from Ti. to Q. 
Then 

limsup Tr A (V^a+|A|g(a a ,b a )\ < gup £)) + 8 (u) - w(Z>)) (6.8) 

In particular u}(G(A, B)) defined as (|2~T7|) exists. 

To appreciate the similarity between (|6.8p and (|3.6J) , one should realize that D 
is a local energy operator, as in (|3.6|) . The proof of this lemma in the case 
that A = B is in the original paper [16 . The proof for the more general case is 
presented in Section [71 Of course, one can prove that the RHS of (|6.8j) is also a 
lower bound: it suffices to copy Section [5l 

By the St0rmer theorem, see (|2.14p , each symmetric state u) on IX is a convex com- 
bination of product states, and since all terms on the RHS of (|6.8[) are affine func- 
tions of u>, it follows that the sup can be restricted to product states. Since, more- 
over, all product states are ergodic, we can replace cu(G(A, B)) by g(oj(A),uj(A)). 
Hence, Lemma 16 . 1 1 implies that 

limsup r-i-r Tr A L-Da+WG(A a ,b a )\ < gup (g(u(A)MA)) + a(w) - w(D)) 

A/°Z d W ^ ' wprod. 

(6.9) 

6.2.1. From the extended PRV to the upper bound. Next, we use (16.91) to formulate 
an upper bound on the quantity 

* Tr A (e-<+l A l G (*- ? A V )) (6.10) 

for A a multiple of V. This means that we have to recall that the lattice sites in 
(|6.9p are in fact blocks. We write A* := A/V and choose 



D 


:= H v 


A 


:= X v 


B 


:= Y v . 
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Then, by the extended PRV, 

§M < sup ( g ^A)MB)) + T^-X{u)-^-MD) 

wprod. on B(A* ) \ I V \ \ V \ 

= sup (g(Lj v (X v ) 7 iv v (Y v )) + -^S{Lj v )--^-uj v (H Vj 

u F onBv \ \v\ \V\ 

where s* indicates that this is the entropy density on the block lattice A*, hence 
it should be divided by \V\ to obtain the density on A. Now, let uj be the infinite- 
volume state obtained by taking a block-product over states toy and let u) be its 
'translation-average', as in Section [U By the conclusions of Section 01 it follows 
that Oj is ergodic and s(uj) = S(luv)- Also, we see that 

A8V\ 



\uj v (X v ) - Q(X)\ < ||X|||SuppX 
\uv(H v ) - Q(E*)\ < r($) 



y\ 

dV\ 



\V\ 

and analogously for Yy . Consequently, we obtain 

®M< sup (g(w(X),u(Y)) + a(w)-w(E*))+0^), V / 7L d 

which proves the upper bound for Theorem l3.2[ since the 0(^)-term can be made 
arbitrarily small by increasing V . 

7. Proof of Lemma HTT1 
Let the state /ia on B A be given by 

with 

Z% (D) := Tr A L-Da+WG^b^ 

Naturally, fi A is the finite-volume Gibbs state that saturates the variational prin- 
ciple, i.e. 

^-\ogZ<i{D) = sup (uj a (G(A a ,B a )) + -±-S(cj a )-w a (D) 

\ A \ u> A onB A V l A l 

= fx A (G(A A ,B A )) + ^-S(n A )-fi A (D) (7.1) 

Our strategy is to attain the 'entropy' and 'energy' of the state fi A via ergodic 
states. For definiteness, we assume that G is of the form 

G(A A , B A ) := [A A } k [B A ] 1 for some integers k, I, 

(which, strictly speaking, is not allowed since G(A Al B A ) has to be a self-adjoint 
operator, but this does not matter for the argument in this section). The general 
case follows by the same argument. 

Since fj, A is a symmetric state on B A (hence translation invariant in the restricted 
sense of Section 0]) , we can apply the construction in Section 0J which yields us 
infinite-volume states ft, and p,. Since we will repeat the construction for different 
A, we indicate the A-dependence in /i^ A ^ and fi^', but remembering that these are 
states on the infinite lattice. They satisfy 

s(/2iA}) = W\ SM (7 ' 2) 
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We have also established in Section U] that is ergodic and that the states 
and p^ A ^ approximate /i A for observables which are empirical averages. However, 
we cannot conclude yet that they have comparable values for G(A,B), except in 
the case where G is linear. Essentially, such a comparison is achieved next by using 
the fact that fj,\ is symmetric. 

Choose a sequence of volumes A n such that along that sequence the RHS of (|7.1[) 
converges. Be assume that p An has a weak-limit, as n / oo, which can always be 
achieved (by the weak compactness) by restricting to a subsequence, of A n . We 
call this limit fi. By construction, it is a symmetric state. 

Energy estimate: Since p An — * p, weakly, and jl Arl (D) = p, An (D), we have 

fx An (D) ^ fi(D) (7.3) 



G-estimates: Using the symmetry of the state /i A , we estimate 

MA (G(A A , B A )) = M® k A ® l B) + ofcA max (||A||, ||B||) fc+i (7.4) 

where the tensor products 

® k A® 1 B := A® . ® A ® B ® . ® (7.5) 

k copies I copies 

denote that all one-site operators are placed on different sites. Since p\ is sym- 
metric, we need not specify on which sites. The error term comes from those terms 
in the expansion of the monomial where two one-site operators hit the same site. 
Since /i is symmetric, we obtain analogously that 

H(G(A,B)) = fi(® k A® l B) (7.6) 

In particular, the LHS is well-defined. Hence, by combining (|7.4p and (|7.6[) , we 
obtain 

^A„(G(iA„,B A J) -» f,(G(A,B)). (7.7) 

For a more general noncommutative polynomial G as defined in Section 12.21 (not 
necessarily a monomial), the convergence (|7.7p follows easily since G(A\ n , B^ n ) can 
be approximated in operator norm by polynomials. 

Entropy estimates: As established in Section IH we have 

pgSfaA) = s(/2 {A} ), for all A (7.8) 

By the lower semi-continuity of the infinite-volume entropy and the convergence 
jl Arl — > /i, we get that 

limsups(/2 {A " } ) < s(jjl) (7.9) 

Hence 

hm -^5(maJ < s(M) (7.10) 

By combining the convergence results (|7.31 17.71 IT. 10[) . we have proven that there 
is a symmetric state fi such that the RHS of (|6.8p with u) = [i is larger than a given 
limit point of the RHS of (|7.1[) . Since the construction can be repeated for any 
limit point, this concludes the proof of Lemma 16. II 
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Appendix A. Proof of Lemma [2~^1 

To prove Lemma 12. 2) it is convenient to introduce an extended framework: Let 
iTu be the cyclic GNS-representation associated to the state u>, 9) u the associated 
Hilbert space and ip G f) u the representant of the state ui, i.e. 

u(A) = {^M^)^, A Gil (A.l) 

The set 7r w (U) is a subalgebra of B{S)J). Let Uj, G Z d be the unitary representation 
of the translation group induced on 7r w (it), i.e. 

U j n u {A)U;=Tr u (T j A). (A.2) 

Ergodicity of ui implies (see e.g. the proof of Theorem III. 1.8 in [T^]) that 

JL J2u. »^ P ^ (A3) 
1 1 ieA 

where P^ is the one-dimensional orthogonal projector associated to the vector ip, 
and A f Z d in the sense of Van Hove. Using (|A.3[) and the translation-invariancc 
Ujip = ip, one calculates 

tt (X A ) tt (Y A ) </> = TAP E U j w(X)U j ,- j irOnU- j >rl> 

1 1 3,3' eh 

— > P^tt(X)P^tt(Y)iP = lj(X)lj(Y)iP 

tv/"1 d 

for local observables X, Y G il. Taking the scalar product with ip, we conclude 
that lu(X\Y\) — y uj{X)uj(Y). The same argument works for all polynomials in 
X\ , Y\ , thus proving Lemma 12.21 Finally, we remark that one can also construct 
the operators 1,7 as weak limits of X\,Y\, as A /* % d (these weak limits are 
simply multiples of identity: uj(X)1,uj(Y)1). This is however not necessary for our 
results. 
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